DISTRIBUTION  STATEMENT  A:  Approved  for  public  release;  distribution  is  unlimited. 


Improved  Ecosystem  Predictions  of  the  California  Current  System 
via  Accurate  Light  Calculations 

Curtis  D.  Mobley 
Sequoia  Scientific,  Inc. 

2700  Richards  Road,  Suite  107 
Bellevue,  WA  98005 

phone:  425-641-0944  x  109  fax:  425-643-0595  email:  curtis.mobley@sequoiasci.com 

Award  Number:  N0001409C0044 
http://www.onr.navy.mil/sci_tech/32/322/environmental_optics.asp 


LONG-TERM  GOAL 

The  goal  of  this  effort  is  to  incorporate  extremely  fast  but  accurate  light  calculations  into  coupled 
physical-biological-optical  ocean  ecosystem  models  as  used  for  operational  three-dimensional 
ecosystem  predictions.  Improvements  in  light  calculations  will  lead  to  improvements  in  predictions  of 
chlorophyll  concentrations  and  other  water-quality  parameters  (such  as  visibility),  and  of  upper-ocean 
thermal  structure,  relevant  to  naval  needs. 

OBJECTIVES 

Currently  available  ecosystem  models  often  use  very  sophisticated  treatments  of  the  hydrodynamics 
(e.g.,  primitive  equation  solutions  in  terrain-following  coordinate  systems  to  obtain  the  advection  and 
upper-ocean  thermodynamics  and  mixing),  fairly  sophisticated  biology  (e.g.,  primary  production, 
nutrient  utilization,  and  grazing),  but  use  grossly  oversimplified  treatments  of  the.  The  optics 
component  of  coupled  ecosystem  models  is  sometimes  just  a  single  equation  parameterizing  the  scalar 
irradiance  or  PAR  (photosynthetically  available  radiation)  in  terms  of  the  chlorophyll  concentration 
and  a  few  parameters  such  as  the  solar  zenith  angle.  Such  simple  light  models  often  fail  even  in  Case  1 
waters,  and  they  can  be  wrong  by  orders  of  magnitude  in  Case  2  or  optically  shallow  waters. 

The  objective  of  this  year’s  effort  was  therefore  to  further  improve  the  previously  developed  EcoLight 
radiative  transfer  model  and  incorporate  it  into  three-dimensional,  time-dependent,  coupled  physical- 
biological-optical  ecosystem  models.  Those  coupled  models  will  be  used  to  evaluate  various  strategies 
for  updating  the  in- water  spectral  irradiance  so  as  to  obtain  accurate  ecosystem  predictions  while 
maintaining  acceptably  fast  computation  times  in  fully  3D  simulations.  In  particular,  I  wish  to  address 
questions  about  how  often  in  time  and  space  the  spectral  irradiance  needs  to  be  updated  when 
modeling  complex  and  dynamic  ecosystems. 

APPROACH 

ONR  previously  funded  (under  contract  N00014-08-C-0024)  the  initial  development  of  EcoLight,  a 
very  fast  radiative  transfer  model  designed  for  use  in  coupled  physical-biological-optical  ecosystem 
models.  During  its  development  and  initial  evaluations  EcoLight  was  embedded  in  an  idealized 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
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does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

2Q  JQ  2.  REPORT  TYPE 

4.  TITLE  AND  SUBTITLE 

Improved  Ecosystem  Predictions  of  the  California  Current  System  via 
Accurate  Light  Calculations 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Sequoia  Scientific  Inc, 2700  Richards  Road  Suite  107, Bellevue, WA, 98005 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distrihution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

8 

3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


ecosystem  model  based  on  a  special  version  of  the  ROMS-EcoSim  physical-biological  model  with  a 
6x6  horizontal  grid  and  periodic  lateral  boundary  conditions.  The  periodic-boundary  version  of 
ROMS-EcoSim  was  chosen  for  development  and  initial  evaluation  of  EcoEight  because  of  its 
geometric  simplicity  and  small  run  times  even  for  multi-year  simulations.  Ten-year  simulations  of  an 
idealized  open-ocean,  Case  1  water,  ecosystem  showed  that  EcoEight  can  replace  the  simple  analytic 
light  model  used  in  the  original  EcoSim  code  with  less  than  a  30%  increase  in  total  run  times.  That 
previous  work  is  described  in  Mobley  et  al.  (2009). 

Given  the  success  of  the  initial  idealized  simulations  with  the  ROMS-EcoSim-EcoEight  code,  the  next 
step  was  to  further  improve  EcoEight  and  imbed  it  into  a  3D  ROMS-CoSiNE  physical-biological 
model  with  boundary  conditions  representative  of  tropical  Pacific  mid-ocean  waters.  The  CoSiNE 
(Carbon,  Silicate,  and  Nitrogen  Ecosystem)  model  is  described  in  Chai  et  ah,  (2002,  2003,  and  2007) 
and  Eujii  et  al.  (2007).  It  was  originally  intended  that  further  evaluations  also  would  be  made  with 
EcoEight  embedded  into  the  NCOM-CoSiNE  model  for  application  to  the  California  Current  System 
and  the  Monterey  Bay  area.  However,  the  present  intention  is  to  apply  the  ROMS-CoSiNE-EcoEight 
model  to  simulations  of  the  Gulf  of  Maine  after  completion  of  the  Pacific  Ocean  runs  now  underway. 

EcoEight  takes  the  following  philosophy.  It  is  necessary  to  solve  the  radiative  transfer  equation  (RTE) 
in  order  to  incorporate  the  effects  of  the  surface  boundary  conditions  and  to  account  for  all  inherent 
optical  property  (lOP,  namely  the  absorption,  scatter,  and  backscatter  coefficients)  effects.  However, 
once  an  accurate  value  of  the  scalar  irradiance  E'o(z,k)  has  been  computed  to  some  depth  Zo  deep 
enough  to  be  free  of  surface  boundary  effects,  it  is  not  necessary  to  continue  solving  the  RTE  to  greater 
depths,  which  is  computationally  expensive.  In  many  cases  of  practical  interest,  it  is  possible  to 
extrapolate  the  accurately  computed  upper-water-column  irradiances  to  greater  depths  and  still  obtain 
irradiances  that  are  acceptably  accurate  for  ecosystem  predictions.  Eikewise,  it  is  not  necessary  to 
solve  the  RTE  at  every  wavelength  in  order  to  obtain  acceptably  accurate  irradiances  at  the  needed 
wavelength  resolution.  Omitting  every  other  wavelength,  for  example,  cuts  the  run  time  by  roughly 
one  half  There  are  several  user-selectable  parameters  in  the  EcoEight  code  that  allow  the  user  to 
define  the  level  of  accuracy  of  the  computed  irradiances.  The  trade  off  is  of  course  between  run  time 
and  accuracy. 

WORK  COMPLETED 

The  initial  version  of  EcoEight  used  in  the  ROMS-EcoSim-EcoEight  simulations  seen  in  Mobley  et  al. 
(2009)  was  written  with  the  goal  of  getting  it  working  correctly  within  the  ROMS-EcoSim  package, 
without  much  time  being  spent  on  streamlining  the  code.  In  particular,  the  input  and  output  for  the 
initial  EcoEight  code  was  tied  to  the  EcoSim  code  in  how  information  was  passed  between  EcoSim 
and  EcoEight. 

Given  the  success  of  the  initial  simulations  shown  in  Mobley  et  al.  (2009),  EcoEight  was  completely 
rewritten  from  scratch  in  Eortran  95  to  (1)  bring  it  up  to  the  standards  of  the  ROMS-CoSiNE  code,  (2) 
remove  various  un-needed  calculations  and,  most  importantly,  (3)  make  the  code  completely 
independent  of  the  particular  physical-biological  model  calling  it.  The  re-written  code  is  called 
EcoEight-S  (for  Ecosystem  Light  Subroutine)  version  1.0.  EcoLight-S  is  a  re-entrant  “black-box” 
subroutine  package  that  can  be  called  by  any  ecosystem  (or  other)  model.  To  incorporate  EcoLight-S 
into  a  particular  physical-biological  model,  the  user  needs  only  write  an  interface  subroutine  to  define 
the  inputs  needed  by  EcoLight-S  to  solve  the  RTE.  All  communication  between  EcoLight-S  and  the 
user’s  program  is  via  Eortran  95  modules.  When  calling  EcoLight-S,  the  interface  subroutine  converts 
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the  physical-biological  model  outputs  (component  concentrations  as  a  function  of  depth  at  a  given  time 
and  grid  location;  other  information  such  as  time,  location,  sky  conditions,  bottom  depth  and 
reflectance  in  shallow  waters;  and  wavelengths  and  depths  where  irradiances  are  needed  by  the  calling 
model)  into  lOPs  and  boundary-condition  information  on  the  format  needed  by  EcoLight-S  to  solve  the 
RTE.  After  solving  the  RTE  for  the  given  lOPs  and  boundary  conditions,  EcoEight-S  returns  the 
spectral  scalar  irradiance  and  ancillary  quantities  (e.g.,  remote-sensing  reflectance,  down-  and  up- 
welling  plane  irradiances,  and  zenith  and  nadir  radiances)  at  the  requested  depths  and  wavelengths. 

The  user’s  interface  subroutine  then  reformats  the  EcoEight-S  output  as  needed  by  the  calling  physical- 
biological  model. 

The  EcoEight-S  inputs  and  outputs  are  fixed  and  well  documented  in  a  Elser’s  Guide  and  Technical 
Documentation  (Mobley,  2010).  Thus  an  EcoEight-S  user  does  not  have  to  consider  the  internal 
workings  of  EcoEight-S,  and  the  EcoEight-S  code  can  be  maintained  and  further  developed  without 
consideration  of  what  models  will  call  it.  EcoEight-S  is  thus  an  RTE  solver  only;  it  does  not  include 
lOP  and  other  sub-models,  or  a  graphical  user  interface,  as  does  HydroEight.  The  usage  of  EcoEight-S 
within  a  physical-biological  ecosystem  model  is  illustrated  in  Eig.  1. 

E.  Chai  and  I  reported  on  our  initial  work  at  an  ecosystem  modeling  workshop  held  at  the  NATO 
Undersea  Research  Center  in  Ea  Spezia,  Italy,  December  14-18,  2009.  The  development  of  EcoEight- 
S  1.0  is  complete  and  the  code  has  been  thoroughly  tested  and  evaluated.  The  Users’  Guide  and 
Technical  Documentation  are  available.  The  EcoEight-S  software  will  be  introduced  to  the  user 
community  at  the  Ocean  Optics  XX  conference  at  the  end  of  September  2010. 

RESULTS 

Most  ecosystem  models  (CoSiNE  in  particular)  use  PAR  as  the  measure  of  how  much  light  is  available 
for  photosynthesis.  It  is  possible  to  compute  PAR  to  the  bottom  of  the  euphotic  zone  in  a  fraction  of  a 
second  of  computer  time,  with  errors  of  no  more  than  a  few  percent.  Eigure  2  shows  PAR  profiles 
computed  by  EcoEight-S  with  varying  degrees  of  run-time/ac curacy  optimization.  The  un-optimized 
run  was  at  5  nm  resolution  from  400  to  700  nm,  with  the  RTE  being  solved  to  50  m  at  each 
wavelength.  The  chlorophyll  profile  is  seen  in  the  left  panel  of  the  figure.  The  center  panel  shows  the 
optimizations  and  run  times.  Parameter  nwskip  defines  the  number  of  wavelength  bands  to  skip 
between  solutions  of  the  RTE.  Thus  in  this  simulation,  nwskip  =  3  means  solve  the  RTE  at  400  nm, 
then  skip  405,  410,  and  415,  and  then  solve  the  RTE  at  420  nm,  and  so  on.  Irradiances  at  unsolved 
wavelengths  are  obtained  by  interpolation  between  the  solved  wavelengths.  Parameter  Fo  =  0.1  means, 
for  example,  that  the  RTE  was  solved  to  the  depth  where  the  scalar  irradiance  was  0.1  (10%)  of  the 
surface  value  at  each  wavelength.  The  Fq  =  0.1  depth  varied  from  more  than  30  m  near  500  nm  to  less 
than  5  m  near  700  nm.  Irradiances  at  deeper  depths  are  obtained  by  extrapolation  based  on  the  lOP 
profde  below  the  last  solved  depth,  as  described  in  Mobley  (2010).  The  right  panel  shows  the  percent 
difference  in  the  unoptimized  (the  RTE  was  solved  at  all  wavelengths  and  to  50  m  depth)  PAR  profile 
and  PAR  computed  with  various  levels  of  optimization.  We  see  that  the  run  time  can  be  decreased 
from  4.81  s  to  less  than  0.5  s  while  still  obtaining  PAR  to  within  4%  of  the  unoptimized  value  down  to 
50  m  (for  these  lOPs).  PAR  was  obtained  to  within  10%  at  all  depths  with  a  run  time  of  only  0.09  s. 

Eigure  3  shows  computed  PAR  profiles  from  an  ecosystem  simulation  in  which  EcoEight-S  was 
embedded  in  the  combined  physical-biological  model  described  in  Mobley  et  al.  (2009).  The 
chlorophyll  profde  used  for  this  figure  was  the  profile  predicted  for  local  noon  of  day  90  (March  30)  of 
a  five-year  simulation  of  Case  1  water  off  the  coast  of  New  Jersey.  The  chlorophyll  values  varied  from 
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about  3.5  mg  m'  near  the  surfaee  to  less  than  0.5  near  100  m  depth.  For  the  eorresponding  lOPs  the 
0.1%  PAR  level  (roughly  the  bottom  of  the  euphotie  zone)  is  at  about  35  m.  The  unoptimized  run  took 
almost  17  s  for  the  solution  down  to  100  m  at  5  nm  wavelength  resolution.  The  run  time  deereases  to 
3.80  s  if  only  every  fifth  wavelength  band  is  solved  {nwskip  =  4),  whieh  gives  an  RTF  solution  every 
25  nm  and  interpolation  to  5  nm  resolution  in  between.  For  this  optimization,  the  PAR  profile  is  the 
same  as  the  unoptimized  value  to  within  0.7%  to  the  bottom  of  the  euphotie  zone  (35  m),  and  within 
5%  down  to  100  m.  Further  optimization  reduees  the  run  time  to  only  0.31  s  with  PAR  still  the  same 
as  the  unoptimized  run  to  within  4%  down  the  35  m,  and  to  within  30%  at  100  m. 

Table  1  shows  the  times  for  HydroLight  5.1,  EeoLight  5.1,  and  unoptimized  and  optimized  EeoLight-S 
runs  for  simulations  of  optieally  pure  sea  water  and  very  turbid  ease  2  water  eontaining  phytoplankton, 
extra  CDOM,  and  “brown  earth”  mineral  partieles.  The  great  deerease  in  run  times  between 
HydroEight  5.1  and  EeoEight  5.1  are  due  to  the  differenee  in  solving  the  full  RTE  vs.  the  azimuthally 
averaged  RTE.  The  differenee  in  EeoEight  5.1  vs.  the  unoptimized  EeoEight-S  is  due  to  the  greater 
effieieney  of  the  re-written  eode.  These  runs  show  that  it  is  possible  to  eompute  PAR  to  the  bottom  of 
the  euphotie  zone  in  both  very  elear  and  very  turbid  waters  to  within  a  few  pereent  of  the  eorreet  value 
(taken  here  to  be  the  HydroEight  values)  in  less  than  one  seeond  of  eomputer  time.  The  optimized 
EeoLight-S  runs  are  all  over  1,200  times  faster  than  the  eorresponding  HydroEight  runs. 

IMPACT/APPLICATION 

Predietive  eeosystem  models  are  playing  an  inereasingly  important  role  in  our  understanding  of  the 
oeeans.  Applieations  of  sueh  models  range  from  predietions  of  water  elarity  for  military  purposes  to 
management  of  eoastal  waters  for  fisheries.  The  ineorporation  of  the  EeoEight-S  model  into  eoupled 
eeosystem  models  will  give  improved  aeeuraey  in  the  predietions  of  primary  produetion  and  related 
quantities  made  by  sueh  models.  As  the  eoupled  models  beeome  more  trustworthy  in  their  predietions, 
they  will  beeome  even  more  valuable  as  tools  for  oeean  seienee  and  aquatie  eeosystem  management. 
We  therefore  expeet  that  EeoEight-S  will  find  wide  use  beyond  its  initial  applieations  in  the  work 
deseribed  here. 

It  should  be  noted  that  ehlorophyll-based  analytie  light  models  underestimate  the  sealar  irradianee  in 
simulations  of  optieally  shallow  waters  with  bright  refieetive  bottoms.  In  sueh  waters,  bottom 
refleetanee  signifieantly  inereases  the  in-water  sealar  irradianee  and  proportionately  affeets  biologieal 
produetivity  and  water-eolumn  heating  rates.  The  EeoEight-S  solution  of  the  RTE  to  a  given  optieal 
depth  is  not  dependent  on  whether  the  lOPs  deseribe  Case  1  or  2  water.  Its  fast  run  times  are  therefore 
retained  in  applieations  to  all  water  bodies.  Unlike  simple  analytie  light  models,  EeoEight-S  ean 
aeeount  for  the  effeets  of  shallow  bottoms  and  is  valid  for  Case  2  waters.  EeoEight-S  also  eomputes 
related  quantities  sueh  as  the  nadir-viewing  remote-sensing  refleetanee  eorresponding  to  the  bio- 
optieal  state  of  the  eeosystem.  This  allows  for  validation  of  eeosystem  model  predietions  using 
satellite  oeean  eolor  radiometry,  without  an  intervening  step  to  eonvert  a  satellite-measured  radianee  to 
a  ehlorophyll  eoneentration  via  an  imperfeet  ehlorophyll  algorithm. 

RELATED  PROJECTS 

This  work  is  a  oontinuation  of  the  EeoEight  development  under  previous  eontraet  N00014-08-C-0024. 
The  present  work  is  a  eollaboration  between  myself  and  Eydia  Sundman  of  Sundman  Consulting 
(embedding  EeoEight-S  into  ROMS-CoSiNE)  and  Eei  Chai  (ROMS  and  CoSiNE  issues,  and  running 
of  large-seale  3D  simulations)  at  the  University  of  Maine.  Chai  is  funded  separately  for  his  work. 
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Eric  Rehm,  a  Ph.D.  student  at  the  University  of  Washington,  is  using  EcoLight-S  as  the  radiative 
transfer  eore  of  an  implieit  inversion  algorithm  that  retrieves  lOPs  from  in-water  measurements  of 
downwelling  plane  irradianee  and  upwelling  radiance. 

The  EcoLight-S  eode  is  also  being  ineorporated  into  the  final  version  of  the  speetrum-matehing  and 
look-up-table  software  deseribed  in  the  report  on  eontract  N00014-10-C-0209.  In  that  software, 
EeoEight-S  is  used  to  eompute  extensive  remote-sensing  reflectance  databases  as  needed  for  inversion 
of  hyperspeetral  imagery  to  obtain  lOPs,  bathymetry,  and  bottom  classification. 
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Any  physical-biological  (P-B)  ecosystem  model 

1 

Water-column  consituent  concentrations  as 
predicted  by  the  P-B  model 

External  environmental  information  (time  and 
location,  v/ind  speed,  etc.) 

X 

_ 1 _ 

User-written  P-B-EcoLight-S  interface  sub 

•  defines  all  lOPs  from  P-B  constituent 
concentrations 

•  defines  surface  and  bottom  boundary 
conditions  (wind  speed,  bottom  refl.,  ect) 

•  sets  flags  for  RTE  solution  options 

routine: 

Reformats  EcoLight-S  outputs 
as  needed  by  the  P-B  model 

X  t 

The  EcoLight-S  subroutine: 

Solves  the  RTE  given  the  lOPs  and 
boundary  conditions 

Returns  the  irradiances,  remote¬ 
sensing  reflectance,  etc.  for  use 
by  the  P-B  ecosystem  model 

Fig.  1.  Usage  of  EcoLight-S  within  any  coupled  physical-biological  ecosystem  model.  The  interface 
subroutine  makes  EcoLight-S  independent  of  the  calling  model.  [The  figure  shows  a  flowchart  of 
calculations  in  a  coupled physical-biological-EcoLight-S  ecosystem  model.] 


Eig.  2.  PAR  profiles  computed  for  various  levels  of  run-time  optimization  in  the  call  to  EcoLight-S. 
With  no  optimization,  EcoLight-S  took  4.81  s  to  solve  the  RTE  at  all  wavelengths  to  a  depth  of  50  m. 
PAR  can  be  computed  to  within  a  few  percent  of  the  unoptimized  value  in  less  than  0.5  s  of  run  time. 
[The  left  panel  shows  the  chlorophyll  vs.  depth  profile  used  in  all  simulations.  The  center  panel 
shows  the  PAR  vs.  depth  profiles  for  various  levels  of  optimization.  The  right  panel  shows  the 
percent  differences  between  the  unoptimized  and  optimized  solutions.] 
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nwskip  =  9;  F„  =  0.10 


400  500  600  700 

wavelength  [nm] 


Fig  3.  Depth  vs.  wavelength  contours  of  the  spectral  scalar  irradiance  Eo(z,F) 
for  various  levels  of  optimization.  Contour  values  are  in  W  m'^  nm'\  The  black  bars  at  the  bottom 
of  the  figures  show  the  wavelength  bands  where  the  RTF  was  solved.  The 
black  symbols  in  the  bottom  two  panels  show  the  depths  to  which  the  RTE  was  solved  at  the  solved 
wavelengths.  The  unoptimized  run  (upper  left)  took  16.88  s;  the  most-optimized 
run  (lower  right)  took  only  0.31  s  while  giving  the  same  PAR  profile  to  within  4%  to  the  bottom  of 

the  euphotic  zone  near  35  m. 


1 


Table  1.  Run  time  in  seconds  for  simulations  of  pure  water  and  turbid  Case  2  water.  For  each 
simulation,  the  times  in  the  left  column  are  for  a  Xeon  2.00  GHz  CPU  and  Windows  XP/SP3;  times 
in  the  right  column  are  for  an  Intel  Core  15  2.40  GHz,  32  bit  CPU  and  Windows  7.  Runs  were  400- 
700  nm  by  10  nm.  PAR  values  for  all  simulations  were  the  same  to  within  3%  at  all  depths  for  the 

respective  simulations. 


model 

pure  water  to  400  m 

0.1%  PAR  depth  >  400  in 
Secchi  depth  =151  m 

tiirhid  Case  2  to  40  ni 

0.1%  PAR  depth  =20  in 
Seeehi  depth  =4.9  in 

llydroLight  v.  5.1 

811.4 

229.8 

462.0 

123.2 

EcoLight  V.  5. 1 

13.0 

4.5 

7.9 

3.0 

EcoLight-S,  no  optimization 

9.3 

2.78 

5.33 

1.61 

EcoLight-S,  with  optimization 

0.61 

0.19 

0.30 

0.09 

